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NOMENCLATURE 


Equation  (22) 
function.  Equation  (32) 
function,  Equation  (33) 
cross  sectional  area.  Equation  (1) 

7l/Aq,  Equation  (8) 

variable  bar  wave  speed.  Equation  (4) 

c/cg,  Equation  (8) 

Young's  modulus 
E/Eo»  Equation  (8) 

Heavyside  Step  Function,  Equation  (7) 
constant;  area  parameter,  Equation  (20) 
constant.  Equation  (23) 

Laplace  transform  parameter,  Equation  (21) 
lot- X) constant;  wave  parameter,  Equation  (19) 
time.  Equation  (1) 
t  c0/5o.  Equation  (8) 
arrival  time  of  wave,  Equation  (34) 
constant.  Equation  (35) 
axial  displacement.  Equation  (1) 

,  Equation  (8) 

transform  of  displacement.  Equation  (21) 
transform  of  strain,  Equation  (26) 
transform  of  particle  velocity.  Equation  (27) 
axial  coordinate.  Equation  (1) 
x/xQ,  Equation  (8) 


icl{v1piq  Ql  *0|  >>  £ 


constant;  modulus  parameter.  Equation  (15) 
[e-O/a.  ,  Equation  (23) 

constant  density  parameter.  Equation  (16) 
mass  density,  Equation  (1) 

,  Equation  (8) 
axial  stress,  Equation  (1) 

CT/Oo  ,  Equation  (8) 
constant,  Equation  (6) 
transform  of  stress.  Equation  (21) 
strain,  Equation  (2) 

ax 

The  subscript  "o"  denotes  quantities  evaluated  xD 
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INTRODUCTION 


The  objective  of  the  work  presented  here  is  to  determine  mathe¬ 
matically  the  elastic  structural  response  as  regards  to  stress,  strain, 
particle  velocity  and  displacements  in  long  rod  projectiles  subjected 
to  at  one  end  to  impact  loading.  This  study  is  limited  to  rod  pene- 
trators  with  a  continuously  varying  inhomogeneity. 

The  dynamic  structural  response  of  projectiles  subjected  to 
impact  loading  from  a  target  is  an  extremely  important  factor  in 
projectile  design.  Generally,  at  ballistic  velocities,  the  response 
is  best  described  by  a  model  accounting  for  dynamic  plastic  and, 
possibly,  hydrodynamic  strains  in  addition  to  the  smaller  elastic 
strains.  However,  at  very  low'  impact  velocities  a  reasonable 
approximation  is  to  view  the  projectile  response  as  totally  elastic. 
This  will  be  the  assumption  followed  in  the  present  analysis.  Repre¬ 
sentative  material  properties  -  such  as  elastic  modulus,  mass 
density  and  yield  strength,  may  vary  within  the  projectile.  It  is  the 
purpose  of  this  study  to  present  theoretical  closed  form  solutions  of 
the  governing  differential  equations  which  include  the  effects  of  such 
variations  or  inhomogeneities.  These  solutions,  while  of  limited 
practical  value,  are  principally  intended  for  use  in  evaluating 
approximate  or  numerical  methods  which  are  designed  to  aid  in  the 
solution  of  realistic  problems.  As  will  be  shown,  the  specific  forms 
of  inhomogeneity  used  allow  a  wide  range  of  material  property 
variations. 

The  propagation  of  elastic  strain  waves  in  inhomogeneous 
has  been  extensively  investigated  in  recent  years.*  This  is  largely 
due  to  the  fact  that  it  is  now  well  established  that  nonuniform  thermal 
gradients  (produced,  perhaps,  during  the  manufacturing  process) 
will  in  turn  produce  a  variable  modulus  and  density  in  a  solid  or 


* 

See  References  1  through  15. 
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hollow  barlto  ^  or  in  a  cylinder.®*  6  The  propagation  of  elastic  plastic 
waves  in  thin  inhomogeneous  bars  was  considered  by  Cristescu,  'and 
Perzyna.  *  ^ 

In  the  present  analysis  it  is  assumed  that  the  projectile  is  an 
elastic  solid  bar  of  variable  cross  sectional  area  using  the  classical, 
or  Love,  bar  theory.1®  It  is  also  assumed  that  the  bar  strikes  the 
target  at  zero  degrees  obliquity  and  thus  that  the  stress  pulse  is 
always  normal  to  the  bar  axis.  The  use  of  the  classical  bar  theory 
will  not  produce  serious  error  if  the  rate  of  change  of  the  cross 
sectional  area  is  small  and  the  propagated  waves  have  wavelengths 
large  compared  to  the  dimensions  of  the  cross  section.  ^  The  modulus 
of  elasticity  and  density  will  vary  continuously  along  the  bar  in  the 
form  of  a  power  law  in  the  spacial  coordinate.  Such  representation 
has  proved  valuable  to  earlier  investigators  of  plane  waves  in  inhomo¬ 
geneous  rods,1*  2»  5*  W)  or  shear  waves  in  infinite  media.17'  18 


Earlier  investigators  have  treated  this  problem  in  several  ways 
ranging  from  purely  analytical  to  wholly  numerical.  Closed  form 
solutions  have  been  obtained  by  Reiss,1®  Datta, 14  Sur,15  Cooper, 
Reddy,  23  Nayfeh  and  Nemat-Nasser.  24  Datta  and  Sur  obtained 
solutions  in  which  the  modulus  varies,  respectively,  as  a  linear  or 
exponential  function  of  the  space  variable  and  the  rod  was  finite. 
Reiss1®  has  presented  a  class  of  closed  solutions  to  seven  problems 
having  inhomogeneous  properties  by  reducing  all  governing  equations 
to  several  canonical  forms  having  easily  obtainable  solutions.  No 
specific  examples  are  given.  The  technique  presented  here  and  that 
of  Reiss1®  produce  identical  results  in  some  circumstances.  Cooper^® 
has  used  the  method  of  progressing  wave  expansions  to  obtain  solutions 
in  the  form  of  an  infinite  series.  For  a  specific  power  law  inhomo¬ 
geneity  in  the  modulus  the  series  terminates  in  one  term.  This  single 
solution  is  of  the  form  obtained  in  the  present  report.  The  latter 
comment  applies  as  well  to  another  approach  taken  by  Reddy.  23 
Nayfeh  and  Nemat-Nasser24  have  presented  a  large  class  of  closed 
form  solutions  as  well,  but  only  for  time  harmonic  waves. 


* 

See  References  1  through  24. 
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THEORY 


Governing  Equations 


Consider  an  isotropic  but  continuously  inhomogeneous  projectile 
(or  bar)  of  elastic  material,  and  neglect  completely  the  Poission 
effect,  radial  inertia  and  bending.  Thus  the  equilibrium  equation 
for  the  linear  elastic  projectile  is: 


j_ 

ax  +  a 

where  C 


=  the  only  non- vanishing  stress  component 


(1) 


ft  =  the(v?.riable)cross  section  area 


=  the(variable)mass  density 


the  displacement 

the  spacial  and  temporal  coordinates 


To  this  is  added  the  stress-displacement  relation: 


CT-  E  — =r 


(2) 


where  E  is  the  modulus  and 
component. 


c>X 


is  the  only  non- vanishing  strain 


Eliminating  the  stress  between  Equations  (1)  and  (2)  yields  the 
Navier  type  relation: 

tfu.  r  e'  _  J _ 

ax  L  IT  *  J  Cc;2 


(3) 


where 

and  a 


C"  =  E/f 

C)  denotes  the  spacial  derivative. 


(4) 
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The  problem  under  consideration  is  that  of  an  infinitely  long  elastic 
projectile,  motionless  prior  to  Thus  the  appropriate  initial 

conditions  subsidary  to  Equation  (4)  are  f Figure  D 


on 


a_  ^  =  ° 

\>  4.0  5  XoiXiOO 


(5) 


For  a  step  input  in  stress  at  X0  »  the  boundary  condition  is 

E  =  CTC  HU1)  at  X=X0  (6) 

where  0^  is  a  constant  amplitude  and  Htt)is  the  Heaviside  Step 
Function  given  as 

HU):0  4oP  -t^-O  (?) 

H(.t)=  »  -Vor  t>0 

Transforming  to  dimensionless  variables  via  the  following 
definitions: 

*=*/*<>  a=o/a 

e=  e/ec  u=  u|/qx0 

c  =  c/C  /\  =  A/Ao  <8> 

t  =  tc/x0  3=  S/to 

Co  =  \/I/Ic 
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The  Equations  (1  to  6)  are  transformed  into 


^5-  cr  s  s  ^ 
zx  /\  9t2 


— »  c  BVA 


r  -  e 

yu 

"bV 


bx 


~  c  3xa  c  L  A  E-» 3,1 


C*  = 


(9) 

(10) 
(11) 

(12) 


U=  O-  ^  foP:l  ^Xioo  ,-fc -O  (13) 
ox 


7)U 

'd'X 


=  HU)  at  x -  \ 


(14) 


It  is  convenient  to  note  that  Equation  (11)  is  a  linear  wave  equation 
with  variable  coefficients  and  is  expressed  only  in  terms  of  displace¬ 
ment.  Thus  provided  the  wave  speed,  c,  is  real,  the  equation  is  of 
the  hyperbolic  type  and  admits  of  propagating  wave  solutions. 
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Solution  Technique 


It  is  Equation  (11),  subject  to  Equations  (13)  and  (14)  which  is  to  be 
solved.  In  general,  two  types  of  solutions  are  available:  numerical 
and  analytic.  A  numerical  solution,  such  as  the  finite  difference 
form  of  the  method  of  characteristics,  can  be  obtained  for  virtually 
any  continuous  choice  of  E,  ^  .1?  An  analytic  solution  can  be 
obtained,  however,  only  when  the  choice  of  these  same  quantities 
leads  to  a  Sturm -Louiville  problem.  Thus  choices  for  E,  ^  ,  and 
A  must  render  Equation  (11)  amenable  to  solution  in  terms  of  such 
orthogonal  functions  as  the  Legendre,  Bessel,  trigonometric  or 
hypergeometric  functions.  Even  then  only  one  function,  the  Bessel 
function,  allows  an  analytic  solution  of  Equation  (11)  to  appear  in 
closed  form  in  terms  of  elementary  functions.  It  is  this  type  of 
solution  which  is  demonstrated  here. 

Using  References  1,  2,  5,  10,  17,  18, 
and  ^  are  made 

E  =  X* 

%=  xK 

and 

«■  io  •, 

4.  1 


the  following  choices  for  E 

05) 

06) 

X  (17) 


where  and  r\  are  constants.  As  will  be  shown  such  a  choice 
allows  solutions  of  Equation  (11)  to  be  obtained  in  closed  form. 


See  References  1  through  18. 
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(18) 


Equation  (11)  becomes 


■&X* 


vu 


2 

where  c  now  has  the  special  form 

C?=  x."  '  =  X 


2  .°i-X  ^,6 


(19) 


and  s  is  set  equal  to  d->.  In  Equations  (18)  and  (19),  it  has  been 
assumed  that  the  area  is  represented  as  a  power  law  in  the  form 


/\<X)  =  X  * 


(20) 


which  represents  a  solid  of  revolution  about  the  x-axis  as  shown  in 
Figure  1.  Several  typical  projectile  shapes  contained  in  this  analysis 
are  shown  in  Table  1. 


TABLE  1. 

Typical  Projectiles  Incorporated  in  Equation  20 


Value  of  k 

0 

Type  of  Projectile 

Right  circular  cylinder 

1 

Ogival  (truncated  paraboloid) 

2 

Truncated  cone 

12 


It  is  appropriate  to  note  here  tha+  equations  similar  to  Equation 
(18)  govern  solutions  to  other  types  r»t  elastic  wave  propagation 
germane  to  munitions  problems.  These  are: 

1.  An  axially  impacted  membrane  shell  with  variable  density  and 

2 

modulus. 

2.  Axial  shearing  of  the  internal  diameter  of  a  cylindrical  bore  in 
an  elastic  half  space.  ^ 

1  O 

3.  A  spherical  blast  wave  in  an  inhomogeneous  medium. 

18 

4.  A  cylindrical  blast  wave  in  an  inhomogeneous  medium. 

Problem  1  is  applicable  to  a  cartridge  case  which  may  be  modeled 
during  extraction  as  a  truncated  conical  membrane  shell  with  variable 
properties.  The  case  receives  axial  impulses  from  both  burning  pro¬ 
pellant  and  vibratory  contact  with  the  bolt.  Problem  2  describes, 
approximately,  the  radial  propagation  of  axial  shear  waves  in  a  target 
plate  due  to  projectile  penetration  at  normal  obliquity.  This  is  part 
of  the  far  field  shear  motion  generated  at  the  shear  interface  between 
projectile  and  plate.  Problems  3  and  4  describe  the  far-field  response 
of  an  inhomogeneous  medium  (such  as  the  earth)  to  detonation  of  con¬ 
fined  or  buried  explosive. 

The  solution  of  Equation  (18)  satisfying  the  initial  conditions  ^ 
Equation  (13)^s  obtained  by  taking  the  Laplace  transform  of  Equation 
(18)  and  employing  Equation  (13).  The  Laplace  Transform  of  Equation 
(18)  is 


where 


gjj  .  _  Js*  U  =  O 

A  X*  x  dx 

00 

Uvx,p'>  =  /  exp(.-pt).vx*;t)<)t 


See  References  2  through  18. 
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The  solution  of  Equation  (21)  is  given  in  Reference  19  as 

u  cx,P  > = Jamo  xa(|  *) + w 

+oi»  /3  ^  O 

Here  p  is  the  transform  parameter,  A(p)  and  B(p)  are  (as  yet)  un¬ 
specified  constants,  and  U(x,  p)  is  the  transform  ofu(x,  t).  Capital 
letters  will  henceforth  represent  the  transform  of  the  corresponding 
lower  case  variable.  In(z)  and  Kn(z)  are  modified  Bessel  Functions 
of  the  order  n  and  of  the  first  and  second  kinds,  respectively.  For 
convenience,  (\  and  A  are  defined  as 


a  = 

2 

(23) 

/3  = 

2-S> 

a 

and 

a/p  =  b- 

-  <.wv°0]  /  L2--^  (24l 

or 

n  = 

-  al/3  =  \j 

-  ch.v<=a]/L^-2J 

Following  convention,  that  is  considering  |3  to  be  greater  than  zero, 
we  choose  in  Equation  (22)  due  to  the  unbounded  behavior 

of  X|\(z)as  J  approaches  infinity  for (3>0.  No  attempt  to  establish 
this  assumption  rigorously  from  the  usual  "boundedness  at  infinity" 
arguments  of  References  2,17,18  for  boundedness  on  either  displace¬ 
ment^  or  stress^’ is  implied.  For  the  case  of  p>*  the  boundary 
condition^  then  determines  BqO  as 


* 

See  References  2  through  19. 
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and  the  general  transformed  solution  for  displacement  is  then  obtained 
from  Equations  (22)  and  (25).  Similarly  the  transformed  strain  dU/dx, 
stress  £  and  particle  velocity  V  are,  respectively. 


ax 

and 

v  --  P  U 


(26) 


(27) 


(28) 


In  deriving  Equation  (25)  it  has  been  assumed  that  n  takes  on  the 
negative  sign  in  Equation  (24).  This  choice  will  be  sufficiently  general 
for  what  follows.  The  restriction  of  0>  0  will  be  discussed  subse¬ 
quently. 

1 7  1 R 

As  shown  in  References1'  • 10  when  n  becomes  half  an  odd  integer, 
the  Equations  (22),  (26),  (27)  and  (28)  may  be  inverted  into  closed  form 
solutions.  This  is  easily  noted  from  the  fact  that  when  n  is  half  an 
odd  interger,  i.  e. , 


rv=+a/p=^+i  • 

<—  • 


<2’> 


17 

P.C.  Chou  and  P.  F.  Gordon,  "Radial  Propagation  of  Axial  Shear 
Waves  in  Nonhomogeneous  Elastic  Media",  Journ.  Acoust.  Soc. 
Amer. ,  Vol.  42,  No.  1,  pp  36-41,  1967. 

18 

E.  Sternberg  and  J.  Chakravorty,  "On  the  Propagation  of  Shock  Waves 
in  a  Nonhomogeneous  Elastic  Medium",  J.  A.  M. ,  Vol.  26,  Trans. 
ASME,  Vol.  81,  Series  E,  1959,  pp  528-536. 
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17  18 

the  Bessel  functions  become  elementary  finite  series  of  the  form  ' 


and 


pro 

K.„^«  K„<*) 


any  o 


(30) 


(31) 


which  have  the  immediate  benefit  that  the  infinite  series  for  U,  ]T  ,  V 
become  the  ratio  of  two  finite  polynomial  series  in  .  Such  ex¬ 
pressions  are  readily  invertable. 

It  is  seen  from  the  above  expressions  that  corresponding  to  each 
j  there  are  an  infinite  number  of  0(  ,  X  and  V\  which  satisfy  each  term 
of  Equation  (28).  In  a  problem  in  which  p(  ,  ^  ,  and  ^  art  specified 
it  is  immaterial,  by  Equation  (30),  whether  the  negative  or  positive 
sign  of  Equation  (23)  is  chosen.  However  the  form  of  Equations  (24) 
to  (27)  was  derived  using  the  negative  sign;  thus  the  negative  sign 
should  be  used  in  practice. 

In  order  to  demonstrate  the  range  of  solutions  possible  for  various 
shaped  projectiles  we  consider  three  examples:  cylindrical,  ogival 
and  conical.  From  Table  II  these  correspond  to  the  area  parameter 
k,  having  values  of  0,  1  and  2  respectively.  Equation  (29)  is  used  as  a 
generating  function  for  all  closed  form  solutions  by  allowing  j  to  take 
on  every  integer  and  zero.  Specifically,  select  projectiles  of  variable 


P.  C.  Chou  and  P.  F.  Gordon,  "Radial  Propagation  of  Axial  Shear 
Waves  in  Nonhomogeneous  Elastic  Media",  Jour.  Acoust.  Soc. 
Amer. ,  Vol.  42,  No.  1,  pp  36-41,  1967. 

E.  Sternberg  and  J.  Chakravorty,  "On  the  Propagation  of  Shock 
Waves  in  a  Nonhomogeneous  Elastic  Medium",  J.  A.  M.  Vol.  26, 
Trans.  ASME,  Vol.  81,  Series  E,  1959,  pp  528-536. 
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elastic  modulus  ,  but  constant  density  A  '-0  • 

functions,  Equation  (29),  have  the  form: 

The  generating 

cylinder 

•H 

1 

**■» 

N 

II 

(32/ 

ogive 

=  2(l+2j)/(2j-l) 

(33) 

cone 

=  4(1  +  j)/(2j  -  1) 

(34) 

TABLE  II 

Closed  Form  Solutions  for  Various  Shaped 
Projectiles:  E  =  E0X^ 


Cylindrical 
k  *  0 

Ogival 
k  =  1 

Conical 

k  =  2 

j _ 

j _ 01 

1 

0  0* 

0  -2 

0  -4 

1  4 

1  6 

1  8 

2  8/3 

2  10/3 

2  4 

3  12/5 

•  • 

3  14/5 

•  • 

3  16/5 

•  t 

j  •  • 

•  • 

00  2 

•  • 

t  • 

CO  2 

•  • 

•  • 

00  2 

-1  4/3 

-1  2/3 

-1  ub 

-2  8/5 

-2  6/5 

-2  4/5 

-3  12/7 

•  • 

-3  10/7 

•  • 

-3  8/7 

•  • 

*  *  j 

•  • 

-00  2 

•  • 

•  • 

-CO  2 

•  • 

•  • 

-CO  2 

a  Homogeneous  cylinder 
k  Homogeneous  cone 
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A  few  values  of  dl  from  these  generating  functions  are  given  in 
Tabic  II.  It  should  be  noted  that  the  usual  homogeneous  solution  for 
the  cylinder  and  cone  are  special  cases.  An  interesting  fact,  which 
may  be  deduced  from  Equations (32)  to  (34)  or  Table  II  ,  is  that  in  the 
limit  as  j  goes  to  infinity,  negatively  or  positively,  the  limiting  closed 
form  solution  is  (*  =  2  or  =  0  .  Note  that  for  /9  =  O  Equation  (22) 
is  singular. 

Two  detailed  examples  of  the  present  solution  technique  for  the 
case  of  a  paraboloidal  projectile  with  a  variable  modulus  are  now  pre¬ 
sented.  As  stated  earlier,  the  choice  of  k  =  1,  \  =  0,  but  Oi  variable 
leads  to  the  formulation  of  Reference  17.  Two  closed  form  solutions 
not  presented  in  that  reference  follow.  These  correspond  to  the  cases 
of  C&  =  6/5  and  -2.  Following  exactly  the  procedure  of  Reference 
17,  using  Equations  (22)  to  (31),  wo  obtain  the  solutions  in  transform 
space  for  the  displacement  as 


for 


o(=  Vo/5  Cn  =  3/2  T  /3 -  2  /  5  ) 

\Kx,p)=  '  x  ^ V]  >  [^  V  *  l)s) 


for 


a-  -a  t/3 --2,n=-'/0 


TT,  «  -s -Wtt  „ 
TKx,p9=p.x  <  ex 


(36) 


P.  C.  Chou  and  P.  F.  Gordon,  "Radial  Propagation  of  Axial  Shear 
Waves  in  Nonhomogeneous  Elastic  Media",  Journ.  Acoust.  Soc. 
Amer. ,  Vol.  42,  No.  1,  pp  36-41,  1967. 
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The  inversion  of  Equations  (35)  and  (36)  requires  only  a  standard 
transform-pair  table  and  repeated  application  of  the  convolution 
(Faltung)  theorem.  The  final  results  are: 


ot  =  fa/5 


and 


UU,t)=  -X-"*  HttrVJexpC-ll** 

-a,]).  +■ 

" t  x*2'?  cot?  ■*(>-«•, 3 )  -*-5A>  x'2/s  1 

XW5  Htt-oi  2xp(r|  Lvv 

-  A.])  -  ([ 

+  Cxws-  x''YY^U-<v>Y 


(37) 


where 


cuxv-IO^O  ?  *=*5 


(38) 

cr  * 
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and 


(39) 


a*00=  U-0/2. 

The  stresses  in  Equations  (37)  and  (38)  and  the  displacement  in 
Equation  (37)  are  displayed  graphically  in  Figures  2  to  4. 

As  can  be  seen  from  Figures  2  to  4  for  the  case  of  =  6/5  the 
displacement  and  stress  are  well  behaved  and  remain  finite  for  all 
values  of  time.  For  the  somewhat  degenerate  case  of  ^  =  -2  the 
stress  is  well  behaved,  tending  to  zero  as  x  increases.  For  the 
same  problem,  however,  we  note  that  the  displacement  at  a  fixed 
spacial  location  goes  unbounded  with  increasing  time.  This  latter 
case  satisfies  the  condition  of  (3>  0  and  the  regularity  of  stress 
given  by  Equation  (40)  which  will  be  presented  later.  It  is  concluded, 
therefore,  that  even  when  the  wave  speed  is  bounded  and  the  stress 
well  behaved  no  assurance  that  the  displacement  is  bounded  may  be 
implied. 


Wave  Front  Information 


By  virtue  of  the  application  of  the  discontinuous  (step)  stress  at 
x  =  1„  and  the  hyperbolicity  of  Equation  (18),  it  may  be  expected  to 
propagate  waves  in  which  the  stress,  strain  and  particle  velocity  are 
discontinuous,  while  the  displacement  remains  continuous.  Such 
elastic  discontinuities  are  well  understood.  Note^  that  the 

"first  characteristic"  or  locus  in  x  - 1  space  of  the  leading  wave 
front  is  given  by 


$ 

See  Reference  20  to  22. 
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Figure  2.  Stress  as  a  Function  of  Time  at  Several 
Locations  in  an  Ogival  Projectile: 

OC  =  6/5 
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TIME 


Figure  3.  Displacement  as  a  Function  of  Time  at  Several 
Locations  in  an  Ogival  Projectile:  =  6/5 
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TIME 


o 


00 


Figure  4 


CO  Tfr 

•  * 

^‘ssaais 


Stress  as  a  Function  of  Time  at  Several  Locations 
in  an  Ogival  Projectile:  =  -2 
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where  physically,  is  the  arrival  time  of  the  leading  wave  at  the 

location  ^  .  The  graph  of  Equation  (40)  is  shown  in  Figure  5.  As 
X  is  allowed  to  approach  infinity  in  Equation  (40)  then  the  time  of 
arrival  at  j^mpis  given  by 


00 


(°°  5  Si2 

(41) 

-  a/u-s)  ■>  $>2 


If  *0^2,  or  the  wave  speed  increases  without  bound,  and  the 

wave  reaches  the  end  of  the  infinite  rod  in  the  finite  time  as  indicated 
by  Equation  (41).  Such  a  situation  will,  for  the  present  report,  be 
disregarded  on  obvious  physical  grounds.  Thus,  as  stated  earlier, 
only  the  case  of  s  ^  2  (or  (3)0)  will  be  treated.  It  should  be  noticed, 
however,  that  the  case  discarded  here,  0,  has  been  successfully 
treated  in  earlier  works.  It  is  suspected,  but  unproven,  that 

the  solutions  obtained  by  inverting  Equations  (22)  to  (28)  are  mathe¬ 
matically  correct  for  certain  time  intervals  in  the  case  of  |9<0.  This 
opinion  is  supported  in  part  by  References  17,  18  and,  apparently, 
disagreed  with  in  Reference  2. 


In  the  case  of  discontinuous  boundary  data,  discontinuities  or 
"jumps"  in  the  first  derivatives  of  displacement  across  the  leading 
wave  may  be  expressed  analytically.  ^2  If  it  is  assumed  that  the  pro¬ 
jectile  is  initially  unstressed  and  motionless  ahead  of  the  leading  wave 
then  the  values  of  stress,  strain  and  particle  velocity  immediately  be¬ 
hind  the  leading  wave  are  given  by 


<7= 

>  -  -  — 

X 

(42) 

and 

3U  _ 

* 

(43) 

» 

3*? 

ii 

X 

(44) 

See  References  2  through  22. 
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1 


Figure  5 


i  •  i  ■>. 

2  3  a  5 

POSITION, X 


Time  of  Arrival,  t,  as  a  Function  of  x  for 
Various  Values  of  S.  (Equation  40) 
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provided  the  boundary  condition  in  the  Equation  (14)  is  invoked.  Tlie 
expressions  in  the  Equations  (42)  to  (44)  apply  only  at  the  time  ta,  in 
Equation  (40),  and  only  for  |3  0.  * 

An  application  of  the  preceding  wavefront  information  theory  is 
given  in  Table  III.  This  shows  the  instantaneous  values  of  stress, 
strain  and  strain  energy  for  the  elastic  projectiles  of  Table  II  at  the 
wavefront.  These  projectiles  have  a  variable  modulus,  E  =  EoX  ^  , 

and  constant  density,  X-0.  (The  coordinate  x  is  the  location  of  the 
wavefront  corresponding  to  the  arrival  time,  ta,  of  Equation  (40)  ).  In 
a  homogeneous  structure  the  stress,  strain  and  strain  energy  are 
either  constant  or  are  decreasing  as  the  wave  proceeds  down  the  pro¬ 
jectile.  This  effect,  constancy  or  decrease  of  the  wave  strength,  is  a 
fcrm  of  geometric  (areal)  attenuation.  As  the  c rossectional  area  is 
either  constant  (cylinder)  or  increases  with  x  (ogive,  cone)  the  stress, 
strain  and  strain  energy  remain  constant  (cylinder)  or  decrease 
(ogive,  cone).  However,  in  the  same  structures  with  inhomogeneity 
the  stress,  strain  and  strain  energy  may  be  markedly  different  as 
shown  in  Table  III.  This  could  be  important  in  estimating  failure. 
Some  modern  failure  theories  postulate  that  failure  is  related  to  strain 
energy.  Failure  in  these  three  homogeneous,  infinite  projectiles  (if 
it  can  occur)  will  occur  at  the  impact  surface  (x  =  1),  which  is  the 
position  of  maximum  strain  energy.  For  an  inhomogeneous  structure 
the  strain  energy  may  actually  be  greater  at  some  point  other  than 
the  impact  location.  Thus  it  cannot  be  assumed  that  failure  can  for 
such  structures  need  be  localized  only  to  the  impact  surface. 

Information  of  the  type  found  in  expressions  in  the  Equations  (42) 
to  (44)  may  be  found  for  other  classes  of  inhomogeneity  than  discussed 


The  question  of  the  applicability  of  Equations  (42)  to  (44)  for  <•  0 
or  *72 is  unresolved  in  the  present  report.  It  was  found  in  Reference 
17  that  unique  solutions  for  S >2.  were  obtainable  by  the  method  of 
characteristics  only  for  certain  restricted  regions  in  x  -  t  space. 
The  Equations  (42)  to  (44)  arise  from  considerations  in  the  method  of 
characteristics. 
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TABLF  Iff 


Distribution  of  Stress,  Strain  and  Strain  Energy 
at  the  Wavefront  for  Various  Projectiles 


INHOMOGENEOUS 

j  Projectile 

iij i n  v.. 

Ogival 

Conical 

Stress 

01/4 

X 

N(ct-  2)/4 

X(C<-  4)/4 

Strain 

x-3<*  /  ) 

x-(30t+  2)/4 

X-(3C(  +  4)/4 

Strain  Energy 

x-  d/2 

x-(a+  21/2 

. 

x-(<*  +  4 )  /  2 

HOMOGENEOUS 

Projectile 

Cylindrical 

Ogival 

Conical 

Stress 

1 

x-1/2 

x-1 

Strain 

1 

x-!/2 

X-1 

Strain  Energy 

1 

x-1 

-2 

X 
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here.  For  the  general  formulation  the  reader  is  referred  to  Refer¬ 
ences  (20)  to  (22). 


As  can  be  seen  from  Eq  lation  (42)  unless 


\*oc<  ah 


(45) 


is  satisfied,  the  usual  regularity  condition,  i.  e. , 

=  O,  ^0  (46) 

cannot  be  satisfied  along  the  wave  front.  Previous,  but  more  restric¬ 
tive,  investigations*' *  have  shown  that  the  type  of  condition  deduced 
in  Equation  (45)  is  equivalent  to  choosing  In  the  present 

paper,  however,  by  recalling  Equation  (23)  one  sees  that  |3>Odoes 
not  necessarily  imply  Equation  (45).  It  would  appear  that  even  iorf&O 
the  regularity  condition.  Equation  (46),  will  not  necessarily  be 
satisfied. 


SUMMARY 


The  principal  result  of  this  study  was  the  successful  application  of 
two  mathematical  techniques  to  the  solution  of  the  problem  of  the 
linear  elastic  wave  propagation  in  an  inhomogeneous  bar.  The  first 
technique,  the  Laplace  transform,  provided  closed  form  expressions 
for  the  states  of  dynamic  stress,  strain  and  displacement  in  the  bar. 
The  modulus,  density  and  cros sectional  area  were  simultaneously 
varied  with  location  in  the  rod.  The  second  technique,  based  on  the 
theory  of  discontinuities,  provided  simple  expressions  for  the  arrival 
time,  dynamic  stress,  strain  and  strain  energy  at  the  wave  front  in 
the  projectile.  The  values  of  dynamic  stress  and  strain  obtained 
from  both  methods  were  found  to  be  identical. 


jjj 

See  References  17  through  22. 
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The  first  technique  v'as  applied  to  a  paraboloidal  projectile  with 
constant  density  and  two  types  of  variable  modulus,  one  increasing 
and  one  decreasing  with  distance  from  the  point  of  impact.  The  results 
for  an  increasing  modulus  projectile  (C{-  6/5)  in  Figures  2  and  3 
showed  that  the  stress  and  displacement  decreased  monotonically  with 
distance  from  the  impact  end  for  all  times.  The  peak  values  of  stress 
at  every  location  occurred  at  the  time  of  wave  arrival,  while  the  peak 
displacements  occurred  a  short  time  after  wave  arrival.  It  was  found 
that  while  the  displacement  was  continuous  for  all  x  and  t,  the  stress 
was  singularly  discontinuous  at  the  wavefront.  The  results  for  a  de¬ 
creasing  modulus  projectile  (  0(  -  -2)  in  Figure  4  and  Equation  (38) 
showed  again  that  stress  and  displacement  monotonically  decrease 
with  distance.  The  stress,  however,  jumped  discontinuously  at  the 
wavefront  from  zero  to  a  steady  ‘bate  value  and  remains  constant 
thereafter.  This  lack  of  transient  behavior  was  completely  different 
from  the  case  of  Oi  =  6/5,  in  which  the  steady  state  values  of  stress 
and  displacement  are  preceded  by  a  transient  behavior. 

In  addition,  the  displacement  for  (X  =  -2  has  no  finite  steady  state 
value,  it  increases  indefinitely  with  time  at  every  location.  It  was 
also  found  that  the  peak  dynamic  stress  is  greatest  in  the  (X  =  6/5 
projectile,  but  that  the  steady  state  stress  was  greatest  in  the  0(  =  -2 
projectile.  At  a  location  .  3  (nondimensional)  units  from  the  impact 
end  (x  =  1.  3)  the  peak  values  of  stress  for  0(  =  6/5  and  0(=  -2,  respec- 
tiv  ly,  are  85%  and  76%  of  the  applied  impact  stress.  The  corre¬ 
sponding  values  of  the  final  steady  state,  or  static,  stress  are,  re¬ 
spectively  67%  and  76%. 

Results  of  the  application  of  discontinuity  theory  (Table  III)  showed 
clearly  that  for  a  homogeneous  projectile  the  stress,  strain  and  strain- 
energy  at  the  wavefront  are  either  constant  (right  cylinder)  or  decrease 
(ogive,  cone)  as  the  wave  proceeds  down  the  rod.  These  results  were 
due  to  areal  attenuation,  i.  e. ,  the  constancy  (cylinder)  or  increase 
(ogive,  cone)  of  area  over  which  the  stress  acted.  For  the  same,  but 
inhomogeneous,  projectiles  the  results  are  markedly  different.  An 
increase  or  decrease  in  stress,  strain  and  energy  at  the  wavefront  de¬ 
pended  upon  the  amount  of  inhomogeneity  (value  of 01  ).  For  example, 
in  an  ogival  projectile  the  stress,  strain  or  strain  energy  increased 
at  the  wavefront  if  the  values  of  were  greater  than  2,  less  than  -2/3 
or  less  than  -2,  respectively. 
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Based  on  these  results,  the  following  conclusions  are  drawn: 


1.  The  characteristics  of  linear  elastic  wave  propagation  in  in¬ 
homogeneous  long  rod  penetrators  are  markedly  different  from  those 
found  in  homogeneous  penetrators.  Specifically,  it  is  found  that  the 
wave  speed  and  arrival  time  are  functions  of  position  and  that  the 
states  of  stress  and  strain  depend  upon  the  amount  of  ii.  homogeneity. 

2.  If  material  failure  is  assumed  to  depend  only  on  stri  ‘.n  energy, 
then  failure  may  occur  at  locations  other  than  the  impact  surface. 

3.  Certain  classes  of  inhomogeneity  led  to  wave  speeds  such  that 
the  wave  front  reaches  infinity  in  a  finite  time.  These  cases,  while 
mathematically  correct,  are  physically  impossible  and  represent 
restrictions  on  the  applicability  of  the  theory. 
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